Multiscaling in Models of Magnetohydrodynamic Turbulence 
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From a direct numerical simulation of the MHD equations we show, for the first time, that 
velocity and magnetic-field structure functions exhibit multiscaling, extended self similarity (ESS), 
and generalized extended self similarity (GESS). We also propose a new shell model for homogeneous 
and isotropic MHD turbulence, which preserves all the invariants of ideal MHD, reduces to a well- 
known shell model for fluid turbulence for zero magnetic field, has no adjustable parameters apart 
from Reynolds numbers, and exhibits the same multiscaling, ESS, and GESS as the MHD equations. 
We also study dissipation-range asymptotics and the inertial- to dissipation-range crossover. 
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The extension of Kolmogorov's work (K41) jO on fluid 
turbulence to magnetohydrodynamic (MHD) turbulence 
yields || simple scaling for velocity v and magnetic- 
field b structure functions, for distances r in the inertial 
range between the forcing scale L and the dissipation 
scale 77^. Many studies have shown that there are multi- 
scaling corrections to K41 scaling in fluid turbulence [ 1 . 
Solar- wind data Q and recent shell-model studies [^](| 
for MHD turbulence yield similar multiscaling. We eluci- 
date this for homogeneous, isotropic MHD turbulence, in 
the absence of a mean magnetic field, by presenting the 
first evidence for such multiscaling in a pseudospectral 
study of the MHD equations in three dimensions (hence- 
forth 3<iMHD). We also propose a new shell model with 
no adjustable parameters (apart from Reynolds num- 
bers) which displays this multiscaling and reduces to the 
Gledzer-Ohkitani-Yamada (GOY) shell model @J| for 
3d fluid turbulence when b = 0. To extract multiscaling 
exponents we develop the ideas of extended self similar- 
ity (ESS ) pd]Jl^| and generalised extended self similarity 
(GESS) |l2| , |l3[| in both real and wave-vector (henceforth 
k) spaces, that have been used in fluid turbulence ||TT[| - 
& 

We use the structure functions Sp = (|a(x + r) — 
a(x)| p ), where a can be v, b, or one of the Elsasser 
variables Z ± = v ± b, x and r are spatial coordinates, 
and the angular brackets denote an average in the sta- 
tistical steady state. Sp ~ at high fluid and mag- 
netic Reynolds numbers Re and Re^, respectively, and 
for the inertial range 20^ < r « L. The extension 
of K41 to homogeneous, isotropic MHD turbulence 
with no mean magnetic field yields — p/3. Shell 
models [|||] and solar-wind data [Q have obtained mul- 
tiscaling in MHD turbulence, i.e., £™ = p/3 — 5Q, with 
8C,p > and Cp nonlinear, monotonically increasing func- 
tions of p. Work on fluid turbulence suggests || an 
extension of the apparent inertial range if we use ESS 
and GESS @): Thus with ESS, in which fol- 



lows from Sp 



[<Sf ]^p'^ , we should expect by anal- 
ogy that it extends down to r ~ 5r?d (as exploited in 
some MHD shell models §,§). In GESS, which em- 
ploys Qp(r) = Sp(r)/[S$ (r)} p / 3 and postulates a form 

W ~ [0J(r)]'5., withjC = [C p a -pQ/3}/[Q -q$/3], 
it has been suggested O] for fluid turbulence that the 
apparent inertial range is extended to the lowest resolv- 
able r; however, fc-space GESS (l2) shows a crossover 
from inertial- to dissipation-range asymptotic behaviors. 
GESS has not been used in MHD turbulence so far. 

Our studies yield many interesting results: The multi- 
scaling exponents we obtain from 3efiMHD and our shell 
model studies agree (Figs, la and lb) and Cp > Cjf > 
Cp > Cp- Cp ne close to the She-Leveque (SL) predic- 
tion for fluids (C^ L = p/9 + 2[1 - (2/3) p / 3 ]), but 
Cp lie below it (Fig. lc) [0. These differences be- 
tween velocity and magnetic-field exponents are also mir- 
rored in differences in the probability distribution func- 
tions (Fig. Id) for 5v a (r) — w Q (x + r) — v a (x) and 
Sb a (r) — & Q (x + r) — b a (x). ESS works both with real- 
and fc-space structure functions (Fig. 2). To study the 
latter we postulate fc-space ESS (for real-space structure 
functions we use S and Q and for their k— space analogs 
(not Fourier transforms) S and G): 

S a p = (|a(k)H « A%(S^\ L- 1 « k < 1.5k d , 

S; = (\a(k)\P)-A a Dp (S 3 r a P , 1.5fc d <fc«A, (1) 

where A^ p and A a Dp are, respectively, nonuniversal am- 
plitudes for inertial and dissipation ranges and A -1 the 
(molecular) length at which hydrodynamics breaks down 
(cf. m for fluid turbulence). We find that a a p ± Q . 

In our shell model C, p a — Cp: but our data for 3dMHD 
suggest C' P a = 2(Cp 1 + 3p/2)/ll (i.e., S£(Jfe) ~ fc"(^+ 3 P/ 2 ) 
in the inertial range fllof ) ; the difference arises because of 
phase-space factors |[L2[. 
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FIG. 1. (a)-(c)Inertial-range exponents versus p from typical 3dMHD and shell-model runs (Table 1) and their comparison 
with the SL formula: (a) Q/Q, (b) Cp/Ca, and (c) Q, Q + , and £*~ from SH2 - ( d ) Semilog (base 10) plots of probability 
distributions P(8v a (r)) and P(Sb a (r)), with r in the dissipation range; a Gaussian distribution is shown for comparison. 
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FIG. 2. Log-log plots (base 10) of Si versus S3 showing k space ESS for 3dMHD with (a) a = v and (b) a — b. Insets 
illustrate real-space ESS for 3dMHD and ESS for our shell model; the lines show the inertial-range asymptotes. 



C p a and a° seem universal (the same for all our runs (Ta- 
ble 1)); a" is close to, but systematically less than, p/3. 
The k dependences of follow from that of S$ . We find 

S£« B?fc^3-9/2 ; L-i 4£ js <i.5kd, (2) 
5|w B a D k s " exp(-c a k/k d ), 1.5fc d <fc<A, (3) 



where Bf and Bp are nonuniversal amplitudes (Eq. (2) 
holds |L2| for 3dMHD; for our shell model the factor 
9/2 is absent). Thus all S£ ~ k e p exp(-c a a£fc/fc d ) for 
1.5k d < k < A, with Q a v = a a p 5 a (cf. g| for fluid tur- 
bulence). In Eq.(3) 5 a ,c a , and kd are not universal. 
However, we extract the universal part of the inertial- 
to dissipation-range crossover via our fc-space GESS. We 
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first define G p 



S£/(S£) p / 3 ; log-log plots of G a v ver- 
sus G q yield curves with universal, but different, slopes 
for asymptotes in inertial and dissipation ranges. The 
inertial-range asymptote has a slope p p q (as in real-space 
GESS ) ; the dissipation-range one has a slope u a (p, q) = 
[at — p/3]/[a q — q/3}. These slopes are universal, but 
not the points at which the curves move away from the 
inertial-range asymptote. To obtain a universal crossover 
scaling function (different for each (p, q) pair because 
of multiscaling) we define \og(H pq ) = .D" log(Gp) and 
log(if° ? ) = D a qp log(G^); the scale factors D a pq 
nonuniversal, 
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FIG. 3. Log-log plots (base 10) of H§ § versus Hg 9 and 
(inset) Hg 6 versus Hg 9 , illustrate our GESS showing the 
universal inertial- to dissipation-range crossover; lines denote 
inertial-range asymptotes. 

but plots of log(Hp q ) versus \og(H qp ) collapse onto a uni- 
versal curve within our error bars for all k, Rex, and Re b x 
(Fig. 3). 

The MHD equations are 



dZ ± 



(Z=F.V)Z d 



2r,± 



;_V 2 Z^ - Vp* 



(4) 



where v± = (y v ± v b )/2, v v and v b are, respectively, the 
fluid and magnetic viscosities, p* is the effective pres- 
sure and the density p = 1, f ± = (f ± g)/2, and f 
and g are the forcing terms in the equations for dv/dt 
and dh/dt. We assume incompressibility and use a 
pseudospectral method JO] to solve Eq.(Q) numerically. 
We force the first two fc-shells, use a cubical box with 
side Lb = 27r, periodic boundary conditions, and 64 3 
modes in runs MHD1 and MHD2 and 80 3 modes in run 
MHD3 (Table 1). We include fluid and magnetic hyper- 
viscosities (i.e., the term — (v v + v V Hk 2 )k 2 in the equa- 



tion for dv(k)/dt and the term —{v b + v b nk 2 )k 2 in the 
equation for dh(]i)/dt, where H stands for hyperviscos- 
ity). For time integration we use an Adams-Bashforth 
scheme (step-size St). We use Rex = v rms X/v v , Re b x — 
b rms X/u b , X v = [f o °° E v (k)dk/ r o °° k 2 E v {k)dk] 1 ' 2 , X b = 
[f o °° E b (k)dk/ r o °° k 2 E b (k)dkYI 2 , E v (k) ~ S%(k)k 2 and 
E b (k) ~ S 2 (k)k 2 . Parameters for runs MHD1-3 are given 
in Table 1, where r ea = LB/a rms is the box-size eddy- 
turnover time for field a and ta the averaging time; initial 
transients are allowed to decay over a period r t . We use 
quadruple-precision arithmetic; results from our 64 3 and 
80 3 runs are not significantly different. 

Shell models for MHD turbulence have been proposed 
earlier |5^,[io|] , but there is no MHD shell model that 
enforces all ideal 3<iMHD invariants and which reduces 
to the GOY shell model for fluid turbulence, when 
magnetic-field terms are supressed. We present such a 
model and show that it yields ( p m agreement with those 
we obtain for 3dMHD. Our shell-model equations 



dt 



ic„ 



k 2 7± 



v -kn z n + fn 



(5) 



use the complex, scalar Elsasser variables z^ = (v n ±b n ), 
and discrete wavevectors k n — k q n , for shells n; 



= [aik n z^ +1 z^ +2 + a 2 k 



a§k n 



(0) 



which ensures z„ , z n 



^,-1/3 jg a stationary solution 
in the inviscid, unforced limit ]8| and preserves the 
v + , Z + «-> z/_, Z~ symmetry of 3<iMHD. We fix five of the 
parameters, a\ — etg, by demanding that our shell- model 
analogs of the total energy (= ^n(b«l 2 + |&n| 2 )/2), the 
cross helicity (= 1/2 J2n( v ^n~^ v n^n)), and the magnetic 
helicity (= Yl,n(~^) n \^n\ 2 /k n ) be conserved if v± = 
and = 0; while enforcing the conservation of energy, 
we also demand Q that the cancellation of terms occurs 
as in 3g?MHD. We fix the last parameter by demanding 
that, if b n — for all n, our model should reduce to 
the GOY model, with the standard choice of parameters 
Q that conserves fluid helicity in the inviscid, unforced 
limit. Thus, apart from the Reynolds numbers, our shell 
model has no adjustable parameters and a\ = 7/12, a 2 = 
5/12, a 3 = -1/12, a 4 = -5/12, a 5 = -7/12, a 6 = 1/12, 
and q = 2. We solve Eq. (^) numerically by an 
Adams-Bashforth scheme (step size St), use 25 shells, 
force the first fc-shell El, set k Q = 2~ 4 , and use E v = 
S^ikn) / k n , \ v = (2tt/ k )[T, n S2(k n )/'S n k 2 l S2(k n )] 1 ^ 2 , 
X b = (27r//c )[S n S , 2(fc n )/E n fc^S , 2(A; I i)] 1 / 2 , v rms = 
[fc E„^(fc„)/^] 1/2 and b rrns = [k j: n S b 2 (k n )/n}^ 2 . Pa- 
rameters for our four runs SH1-SH4 are given in Table 
1. These use double-precision arithmetic, but we have 
checked in representative cases that our results are not 
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affected if we use quadruple-precision arithmetic. As 
in the GOY model the structure functions S p (k n ) oscil- 
late weakly with k n because of an underlying three-cycle 
I&IJ. These oscillations can be removed either (a) by 
using ESS plots or (b) by using the structure functions 
K, P = <3[a«an+ian+2 + a n _ia„a„ + i/4]P/ 3 ) @. Method 
(a) yields the exponent ratios (Cp/Cf)) which we find are 
universal. Method (b) gives exponents Cj- These have 
a mild dependence on Re\ and Reb\ but this goes away 
if we consider the exponent ratios (p/CSi as m the GOY 
model thus the asymptotes in our ESS and GESS 

plots have universal slopes. 

The Navier Stokes equation (3dNS) follows from 
3<iMHD if we set b = or, equivalently, Reb\=0. How- 
ever, if we start with Rebx — 0, the steady state is char- 
acterised by the MHD exponents and Re\/Reb\ — 0(1) 
(i.e., an equipartition regime) jL8|. Since our MHD shell 
model reduces to the GOY model as Reb\ — ► 0, we 
use it to study the fluid turbulence to MHD turbulence 
crossover, instead of doing costly pseudospectral studies: 
A small initial value of Rebx yields a transient during 
which we obtain GOY-model exponents, but eventually 
the system crosses over to the MHD turbulence steady 
state ||. 

In conclusion, then, we have shown that structure func- 
tions in 3eflVlHD turbulence display multiscaling, ESS, 
and GESS, with exponents and probability distributions 
P(5v a (r)) and P(5b a (r)) different from those in fluid tur- 
bulence. Our new shell model (a) gives the same multi- 
scaling exponents as 3dMHD and (b) reduces to the GOY 
shell model as Re b \ -»• 0. Our ESS and GESS stud- 
ies help us to uncover an apparently universal crossover 
from inertial- to dissipation-range asymptotics. It would 
be very interesting to compare our results with experi- 
ments on MHD turbulence, but two points must be borne 
in mind: (1) solar-wind data might yield multiscaling ex- 
ponents different from ours because of the presence of 
a mean magnetic field; (2) the crossover from inertial- 
to dissipation-range asymptotics might not apply to the 
solar wind because a hydrodynamic description might 
break down in the dissipation range |l9| . However, our 
results should apply to MHD systems which show an 
equipartition regime It would also be interesting to 
see whether the agreement of fjj with the SL formula is 
fortuitous or significant. 
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TABLE I. The viscosities and hyperviscosities v v ,v b , v vH and v bH , the Taylor-microscale Reynolds numbers Re\ and Re b \, 
the box-size eddy-turnover times r e „ and r e b, the averaging time ta, the time over which transients are allowed to decay n, 
and kd (dissipation-scale wavenumber) for our 3dMHD runs (k m ax = 32 for MHD1 and MHD2 and k ma x = 40 for MHD3) and 
shell-model runs SH1-4 (k max = 2 25 k ). The step size(<5i) is 0.02 for MHD1-3, 2.10" 5 for SH1-2, and 10~ 4 for SH3-4. Note 
that r ev ~ 8ti the integral time for our MHD runs. 



Run 


v v 


VvH 


v b 


VbH 


Rex 


R&m\ 


T ev /St 


T ev /6t 


T~t/T ev 


ta/t bv 


k>max j 


MHD1 


8 ■ 10" 4 


7- 10" b 


10~ 3 


8- 10" b 


~ 24.8 


~ 14.3 


~ 8.8 • 10 3 


~ 6 • 10 3 


~ 2 


~ 2.3 


~ 1.83 


MHD2 


8 ■ 10~ 4 


9- 10" 6 


8 ■ 10~ 4 


9- 10~ 6 


~ 24.1 


~ 18.1 


~ 8.8 • 10 3 


~ 5.6 • 10 3 


~ 2 


~ 2.3 


~ 1.83 


MHD3 


8 ■ 10~ 4 


9- 10~ 6 


8 ■ 10~ 4 


9- 10~ 6 


~ 26 


~ 19.6 


~ 7.9 • 10 3 


~ 4.8 • 10 3 


~ 1 


~ 2.2 


~ 2.22 


SHI 


io- 9 





10- 9 





~ 4.6 • 10 8 


~ 7.8 • 10 8 


~10 7 


~ 6 • 10 6 


~ 50 


~ 450 


~ 2 5 


SH2 


io- 8 





io- 8 





~ 4.3 • 10 7 


~ 6.5 • 10 7 


~10 7 


~ 6 • 10 6 


~ 50 


~ 450 


~2 8 


SH3 


10~ 6 





2 ■ 10~ 6 





~ 4 • 10 6 


~ 3 • 10 6 


~ 2 • 10 6 


~ 10 6 


~ 500 


~ 2500 


~ 2 10 


SH4 


4 ■ 10~ 6 





10~ 6 





~ 1.2 • 10 5 


~ 1 • 10 6 


~ 10 6 


~ 1.7- 10 6 


~ 500 


~ 3000 


~2 U 
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